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Abstract 



We show that the isentropic subclass of Buchdahl's exact solu- 
tion for a gaseous relativistic star is stable and gravitationally 
bound for all values of the compactness ratio u[= (M/R), where 
M is the total mass and R is the radius of the configuration in 
geometrized units] in the range, < u < 0.20, corresponding to 
the ref/n/ar behaviour of the solution. This result is in agreement 
with the expectation and opposite to the earlier claim found in 
the hterature. 
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1. Introduction 



Inspite of non linear and coupled differential equations, var- 
ious exact solutions of Einstein's field equations for static and 
spfierically symmetric perfect fluid spfiere are available in the lit- 
erature [1]. But the problem of finding solutions corresponding 
to a regular (i.e. the finite positive density at the origin which 
decreases monotonically outwards) and gaseous (i.e. the den- 
sity also vanishes at the surface of the configuration together 
with pressure) sphere is rather more difficult. Such solutions 
are, therefore, severely restricted to a limited number. To the 
knowledge of the present author, there exist only two exact so- 
lutions in this category which can serve as realistic models of 
relativistic stars. The first solution is Tolman's type VII solution 
with vanishing surface density [2, 3] and other one is Buchdahl's 
gaseous model [4, 5, 6]. Since such solutions may also find appli- 
cation to construct the models of relativistic supermassive stars 
[7], together with the usual models of neutron stars [8], their 
investigation is important. Even more important is the study 
of their gravitational binding and stability under small radial 
pulsations. 

The stability of Tolman's type VII solution with vanishing 
surface density has been undertaken in detail by Negi and Dur- 
gapal [7, 9] and they have shown that this solution remains 
gravitationally bound (for all permissible values of u < 0.3861) 
and pulsationally stable (for all values of u < 0.3428). The 
pulsational stabihty of Buchdahl's gaseous solution was studied 
by Knutsen [6] and he found that this solution turns out to be 
unstable for the range of compactness ratio < u < (1/6) (the 
range in which the Buchdahl's gaseous solution satisfies the reg- 
ularity and causahty conditions). However, this result seems to 
be unhkely, because from the study of the behaviour of this so- 
lution, one finds that the solution becomes a classical polytrope 
of index 1 as -u ^ 0, whereas when u = (1/6) the solution is 
equivalent to a sphere of homogeneous energy-density in the core 
(Schwarzschild core) surrounded by an envelope which becomes 
polytropic {n — 1) as the surface is reached. Since both of these 
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solutions, the classical polytrope of index 1 and the homoge- 
neous density sphere, represent the stable configurations in gen- 
eral relativity [10], it looks surprising that Buchdahl's gaseous 
model which becomes one of these solution in the limiting case, 
— > 0, and the 'effective' core-envelope model comprising both 
of these solutions as m (1/6) is unstable for < < (1/6) ! 

In order to resolve the inconsistency between the reasoning 
presented above and the results of Knutscn [6], we also investi- 
gate the pulsational stability of Buchdahl's gaseous model on 
the basis of variational method [11,12] by using the general 
form of the 'trial function', frequently used in the literature 
[10,13] (whereas Knutsen [6] has used only a particular form of 
the 'trial function' which also satisfies the boundary conditions). 
For comparison, we also employ the particular form of the 'trial 
function' used by Knutsen [6]. Furthermore, we also consider 
the gravitational binding of the configuration which is necessary 
for the existence of a realistic structure (Knutsen [6] did not con- 
sider the gravitational binding of the structures). We perform 
our calculations for the range < it < 0.20 (we extend the range 
of our calculations for the interval (1/6) < u < 0.20 where the 

adiabatic sound speed, Vs{= \J\dP/dE)), exceeds the speed of 
light, c. Actually, we are interested to investigate the stability 
for the full range in which the solution behaves like a regular 
solution and not just for the interval in which the 'causality' is 
not violated. Since for the case of high densities, the adiabatic 
sound speed, f^, docs not necessarily represent the 'true' signal 
propagation speed so that the 'causality' may not be violated 
even if > c [14]). Our calculations show that perfect fluid 
subclass of Buchdahl's gaseous solution is stable and gravita- 
tionally bound for all values of u in the range Q <u < 0.20. 

2. Metric and Expressions for Buchdahl's Gaseous 
Model 

The interior geometry of static spherically symmetric per- 
fect fluid configurations in canonical coordinates is given by the 
metric 
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ds"" = e'dt'' - e^dr"" - r'^^dO'' + r''sui^ed<i?) (1) 

where v and A are functions of r alone. Recalling that we 
are using the 'geometrized units' [G = c = 1; where G is the 
Universal gravitation constant and c is the speed of light in 
vacuum], the metric coefficients of Eq.(l) for Buchdahl's gaseous 
solution yield in the following form [5,6] 

e'' ={l-2u)'—^ (2) 



The expressions for pressure, P, and energy-density, now 
read 

SvrPP^ = "'^^""f, (4) 

(1 -2m) (c + a)2 ^ ' 

(1 -2m) (c + a)2 ^ ^ 
where the variables a, 6 and c are defined as 

d = 2u{sm{z) / z) (6) 

b = 2u cos 2; (7) 

c = 2(1-m) (8) 



and 2; is given by 



< Z < TT (9) 



[l + u{sin{z)/z){l - 



^We use notation for parameters, different from those used in [5] and [6], for conve- 
nience 
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here y is the dimensionless radial coordinate (r/i?) (i.e. the 
radial coordinate measured in units of configuration size) varies 
in the range < 7/ < 1. 
Eq.(9) yields 

T^dy —\i-\-ucosz{l — u)~^]dz (10) 
The adiabatic sound speed, Vg^i is given by 

Since the Eqs.(2)-(5) yield the pulsationally stable structures 
(section 4) for all permissible values of one can construct the 
models of physical objects (i.e., neutron stars and supermassive 
stars) on the basis of Buchdahl's gaseous solution. For an as- 
signed value of the compactness ratio, and a suitable value of 
the central energy-density, Eq (the value of at r = 0) , the ra- 
dius i?(and hence the total mass M(= Ru)) of the configuration 
can be obtained by using the following relation 

2 ^ 7r^^(2-5i.)(l-^)^ 

87r£;o(l - 2ii) ^ ' 



4. Dynamical Stability of Buchdahl's Gaseous Model 

The dynamical stability of Buchdahl's gaseous model may 
be assured by using the variational method [11,12] which states 
that the sufficient condition for the dynamical instability of a 
mass is that the right-hand side of the following equation 



^2 f''e^^^-''^/\P + Eyedr^ 
Jo 

Jo 
Jo 
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Jo 

+871 re^^^+''^/^P(P + EyC^dr (13) 
Jo 

vanishes for some chosen 'trial function' ^ which satisfies the 
boundary conditions 

^ = at r = 0, (14) 

and 

= at r = R (15) 

The system, consists of gaseous configurations discussed in 
section 2, is considered fully 'isentropic' so that the quantity 
7 = [{P + E)/P]{dP/dE) = Ti appearing in Eq.(13) repre- 
sents the adiabatic index. The prime appearing in Eqs.(13) and 
(15) denotes radial derivative, 5P is called the 'Lagrangian dis- 
placement in pressure', a is known as the angular frequency 
of pulsation and R is the size of the configuration. In order to 
solve Eq.(13) for the system under consideration, we employ the 
fourth-order Runge-Kutta iteration by using the 'trial function' 
as a power series of the form ^ — bir{l + -\- a2r^ -|- a^r^)e'^/'^, 
(where ai,a2, and as, are arbitrary constants) [10,13]. The pos- 
itive values of the squared angular frequency of pulsation, 
for all the permissible values of u indicate the existence of pul- 
sationally stable configurations in the range < m < 0.20 as 
shown in Table 1. For comparison, we repeat the calculations 
for the choice of the 'trial function' £^ = re^^^ used by Knutsen 
[6]. The positive values of the squared angular frequency of pul- 
sation, (J 2-, for all values of u as shown in Table 1 confirm the 
above finding. We also calculate the fractional binding-energy 
a[= [Mr — M)/M\ where M.^ is called the rest-mass of the en- 
tire configuration ([15])]. The positive values of a for all the 
above-mentioned values of u show the existence of gravitation- 
ally bound structures. 
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6. Results and Conclusions 

The present study shows that the isentropic, perfect fluid 
subclass of Buchdahl's gaseous model is pulsationally stable 
and gravitationally bound for the range of compactness ratio 
in which the solution satisfies the condition of regularity. This 
finding might be helpful for the future investigation of the rela- 
tivistic supermassive stars as models for quasars on the basis of 
this solution, together with the usual models of neutron stars. 
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Table 1: The angular frequency of pulsational, ai, for Buchdahl's gaseous 
models as calculated by the variational method (Eq.l3) for the choice 
of the 'trial function' ^ = 6ir(l + air^ + a2r^ + a3r^)e''/^, where ai = 
(1/I0i?^),a2 = (l/5-R^),a3 = (3/lOi?^), are arbitrary constants, is shown 
in column 4 for different values of the compactness ratio u (column 1) in 
the range < u < 0.20. For comparison, the corresponding angular fre- 
quency of pulsational, f72, calculated for the choice of the 'trial function' 
^ = re"^/^ in Eq.(13) is shown in column 5. The binding-energy per unit 
mass, a[= {Mr — M)/M, where Mr is the rest mass of the configuration] 
is shown in column 3. Column 2 shows different values of the constant 
k[= (1 — 2u)], used by Buchdahl[4] and Knutsen[6] for corresponding val- 
ues of u. 
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